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Abstract. We present an algorithm which computes the Landau constant up to any 
given precision. 



1. Introduction 

For most of the mathematical constants used today, like tt, e or £, very efficient algorithms 
to approximate these constants are known. There are, however, rare cases such as Bloch's, 
Landau's or the Hayman-Wu constant, where conjectures and rough bounds on such con- 
stants are known but the exact value is unknown. Recently, an algorithm to approximate 
Bloch's constant was given in |Rettingcr 2008 1. We continue this line of research and present 
a similar algorithm to approximate Landau's constant up to any precision. 

Because of the close connection between the Landau and Bloch constants we can reuse 
many of the ideas (and motivations) on the computation of the Bloch constant in this 
paper. A core component of the algorithm, namely the computation of maximum discs 
in the image of a (normalised) holomorphic function, is totally different. So injectivity 
does not play any role any longer. What is more, however, homotopy methods used in the 
algorithm for Bloch's constant do not work any longer. This is due to the fact that for discs 
which need not be schlicht, several overlapping parts of the pre- image come into play. To 
decide whether the images do indeed overlap is in general not even computable. Though 
the ideas of the algorithm are inspired by Type-2 theory of effectivity, we formulate our 
results totally independently of this theory. In this way we hope that our result is accessible 
to a wider audience. Finally, we improve even on the parts which could be taken literally 
from the algorithm for Bloch's constant to narrow the gap between the theory and possible 
implementations. 

Landau's constant (see [Landau 1925) ) gives a quantitative version of the fact that non- 
constant holomorphic functions are open. More precisely, it states that for any r > and 
any holomorphic function / defined on a disc IDV(zo) := {z G C | \z — zq\ < r} with f'(zo) / 
there exists a disc of radius |/'(^o)| • r ■ c inside the image f(3 r (zo)), where the constant 
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c > does not depend on /! Obviously c is bounded from above; thus the supremum, the 
so called Landau constant A, exists. The best upper bound known for A, 



is at the same time conjectured to be the exact value of A. However the best lower bound 
known so far is 

, , , 

(see [Robinson 1938] and [Rademacher 1943] ). 
Putting this in decimal representation gives 



i.e. all we know is the constant up to 4 • 10 . 

In this paper we will give an algorithm to compute Landau's constant A up to any 
precision in the sense that on input n G N some rational number q with \q — A| < 2~™ can 
be computed. 

The main idea of our algorithm is to compute for several normalised functions the 
corresponding A-values. Following the definition, it seems that we have to take the infimum 
for all normalised functions, which could not be done in finitely many steps. We will 
overcome this problem by a compactness argument in Section 5. In the next section we 
recall a few notations and fix the machine model which we will use throughout this paper. 
In Section 3 we introduce a subset of all holomorphic functions such that A is already 
determined by this smaller, compact class. In Section 4 we will then show how A/ of a 
single function can be approximated and, finally, in Section 5 we prove the computability 
of Landau's constant. 



We begin this section with a few remarks on the model underlying our algorithms. Our 
main algorithm (Section 5) will only compute on finite words, where we assume some 
straightforward, standard encodings (representations) of the dyadic numbers, i.e. numbers 
of the form m ■ 2 n with m, n G Z, by such words. We do not need any operations on infinite 
words or other structures like R or C, and therefore our algorithms could be implemented 
by classical Turing machines. If some value is given as input to our machine, we also say 
that the machine computes something on this value. In addition, to simplify things, we use 
a second model for intermediate results where we allow infinite sequences as input. One 
can think of this model as a classical Turing machine which can in addition ask some kind 
of oracle for single elements of the input. The output, however, will always be a finite 
word which will be returned after a finite number of steps. This second model is only used 
for simplification reasons, and our final algorithm to compute Landau's constant will not 
depend on this model. 

Let H r (zo) ■= {z G C | \z — z§\ < r} denote the open disc of radius r with centre zq 
and let H r {zo) denote its topological closure, i.e. B r (zo) ■= {z £ C \ \z — zq\ < r}. To 
simplify notation we use D r := D r (0) and D = Di. A normalised (holomorphic) function on 
a domain D, G D, is a holomorphic function / with f'(0) = 1. The space of normalised 



A < 



r(|)-r(f) 




0.5 < A < 0.54325... 
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functions on D is denoted by M(D). For given domain ACCwe denote the radius of the 
largest disc in A by 1(A), i.e. 

1(A) = sup{r | 3z.B r (z) C A}. 

Given a holomorphic function / with D C dom(/) let A/(L>) denote the radius of the largest 
disc in f(D), i.e. 

\ f (D) = l(f(D)). 

Finally let \f := A/(B) and the Landau constant A be the infimum of all \f with / G N(D). 
Obviously, A is already the infimum of all A f for normalized / with the additional condition 

/;<>;• o. 

Throughout this paper we use most of the time only very basic results of complex 
analysis covered by most textbooks (see e.g. [Ahlfors T9 66] ) . The only somewhat more 
advanced result from complex analysis is Koebe's 1/4 Lemma (see e.g. [Rudin 19 87J): 

Theorem 1 (Koebe's 1/4 Lemma). Let / be an injective holomorphic function on some 
disc D r . Then D r .| / , (0) |. 1/4 (/(0)) C /(D r ). 

The representations of the objects we use throughout this paper are introduced next, 
where we implicitly use some kind of efficiently computable pairing function without further 
mention. 

Let Y denote the class of dyadic numbers, i.e. 

Y := {m-2 n \ n,m G Z}. 

We could use rational numbers as well but having implementations in mind, we stick to the 
efficiently implementable dyadic numbers. Identifying complex numbers z = x + i ■ y and 
pairs (x,y), where we denote the imaginary unit by i, we can represent elements of the set 
Y[l] of complex dyadics simply by pairs of dyadic numbers. In this way, it should be clear 
that there are efficient algorithms to approximate operations like •,/,+, — on the complex 
dyadics, in the sense that, given complex dyadics z, z' 7^ and some n G N, we can easily 
compute a dyadic number y such that \y — zjz'\ < 2~ n etc. To simplify notation, we say 
that we can compute a value z up to precision 2~ n if we can, on input n and z, compute a 
(complex) dyadic y such that \y — z\ < 2~ n . If we can compute z up to precision 2~ n for 
all n, then we say that we can compute z up to any precision. 

Let furthermore I denote the class of intervals [c, d] x [d, d'\ where c, d , d, d' are dyadics. 
We consider these intervals as sets of complex numbers. Let II be a finite alphabet. Then 
II* denotes the set of finite sequences (words) over II, e the empty sequence, and IT 00 the 
set of infinite sequences over II. As usually, we think of IT 00 as a topological space where 
{wH°° I w G II*} is a basis of the topology. It is well known that H°° with the above 
topology is a compact space. 

In our algorithm we need to approximate holomorphic functions in a certain function 
space. We will do this by representing this function class by H°° for II = {1,2,3,4}. To 
simplify things we will first define, for given interval I G I, the representation ^ff of the 
complex numbers in /. To this end, let for a given interval I = [c,d] x [d, d'\, the intervals 

[c, c"] X [d, d"] , [c", d] X [d, d"] , [c", d] X [d", d'} , [c, c"] X [d", d'] , 

where c" = (c + c')/2 and d" = (d + d')/2, be denoted by I\, h, h and I4, respectively. 
With these preliminaries, we use sequences ao, or, ... of numbers in II := {1,2,3,4} to 
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represent complex numbers where we define : II* — > I by 

$/H = (*j(w)) a 
for all I G I, u; G II* and a G II, and finally 

*f (a ai . . .) = P| ^/(aoai •••««), 

i 

i.e. we define complex numbers by suitable nested intervals. 

Let now m 1; m,2, . . . denote a sequence of positive dyadics such that m i ' £% con " 

verges for all e G [0, 1). Then the set of sequences 01,02,. . • of complex numbers where 
the absolute values of the real and imaginary parts of a% are bounded by m% for all z > 1 
defines a class F mit m 2j .., of holomorphic functions on D by identifying a sequence a\, 02, • • ■ 
with the power series 1 + Yli>i a i~ z 1 - We are interested only in power series where the first 
coefficient is 1 because this class represents exactly the derivatives of normalized holomor- 
phic functions. We use \I/ to define a representation of these holomorphic functions: Let 
t\ , t% , . . . be a sequence of natural numbers such that any natural number is encountered 
infinitely often, e.g. 0, 0, 1, 0, 1, 2, 0, 1, 2, 3, 0, . . .. Then for a given sequence az\, 02, ■ ■ ■ 
with Qj 6 II and any n G N the sequence t\, £2 ? • • • defines the subsequence a ni , a n2 , . . . of 
all those elements a nj such that t n . = n. In this way we can define a mapping 

.TrOO . TTOO . 771 

T mi,m2,- ' r mi,m,2,... 

by 

** 1)ffl3> .„(«i«2...) = 1 + J^ m „, mn] x[- mn , (% ]K 1 an s ...)^ n 

n>l 

To simplify notations we denote, for given holomorphic function / : D — > C, the anti- 
derivative of /, which maps to 0, by / /, i.e. / / : D ->■ C, (J* /)(0) = and (/ /)' = /. 
Notice that J m2 (0102 • • •) G jV(D) with the above settings. A proof of the following 
theorem can for example be found in [Mueller 1987] or |Retti nger 2008b]. 

Theorem 2. Let a±, 02, ■ ■ ., mi, m%, . . . and ti, t%, . . . be as above. Furthermore let bi, 
62 , ... be a sequence of positive dyadic numbers such that with 

/ ^ mi,rri2,... {^1^2 • ■ •)> 

for all n > 1 we have 

teeBi_ 2 -».|/(;s)| <b n . 

Then, on input r G Y D (0, 1), z G Y[t] and 

ai,a 2 ,..., mi,m 2 ,..., ti, *2, • • - , &i, &2, • • ■ and n G N, 
we can compute J f(z), f(z) and f^ n '(z) on D up to any precision. 
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3. A-bounding Functions 

Similar to the /3-bounding functions in [Rettin gcr 2008 1 we define here a class F\ of holomor- 
phic functions which determines A in the sense that A = inf /eFA A/. Unlike |Rettinger 2008] 
we use different bounds on the coefficients of the corresponding power series and, further- 
more, use the compactness of IT 00 together with the mapping rather than compactness 
of the class F\ itself. 

To start with we will prove suitable upper bounds on (the derivative of) functions / 
whose A- values approximate A: 

Lemma 1. Let c > 1 be given. Then for any / G AA(D) and w G D with |/'(u;)| > 
c/(l — |u>| 2 ) we have 

A/ > c- A. 

Proof. Let w G and / G W(D) with \ f'(w)\ > c/(l - \w\ 2 ) be given. Let g : D -)• D be 
the automorphism 

«4 (2: + w)/(l + to). 
Then / o 5 (0) = /(«,), (/ o 5 )'(0) = f(w) • (1 - \w\ 2 ), and 

f'(w) ■ (f - |u>| 2 ) 

is a normalised function. Therefore 

- • Afo 9 > rjr, r, — jz , — r^r • ^fog > Aft > A. 

c I /'(it;) I • (1 — \w\ z ) 

As /(D) = f o g(D) the statement of the lemma follows. □ 

To simplify things we fix some value c > 1, say c = f + 2~ 100 for the time being. To 
find A it suffices to consider all Xf of all normalised functions / with |/'(^)| < c/(l — \z\ 2 ) 
for all z G D. This immediately gives us a bound on the coefficients of the corresponding 
power series: 

Lemma 2. Let / be a normalised function such that |/'(z)| < c/(l — \z\ 2 ) for all z£l. 
Then f'(z) = Y2 n a n ■ 2 n with 

(a) ao = 1 and 

(b) |a n | < c- e- (n + 2)/2 for all n > 1. 
(Here e denotes the Euler constant.) 

Proof. Item (a) is obvious because / is assumed to be normalised. 
By the Cauchy inequality we have for given n and r G (0; I) 

1 < 1 sup |/'(z)| < 1 
|z|=r 



r n r n 1 — r z 



Choosing r := y/n/(n + 2) we get 

1 c / 2\? n + 2 n + 2 ,_, 

r n 1 — \ n J 2 2 

Let, for n > 1, m n be some dyadic approximation such that 

c • e • (n + 2)/2 < m n < c • e • (n + 2)/2 + 2~ n . 
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Then we have for given e E (0; 1) 

l + E„>l^n-e" < f •(E„e n + En(™+l)-£ n ) + En2-"-C-e 

< ¥-(l^i + (li F ) + (2-c-e). 
A similar bound does also hold for the corresponding derivatives: 

En — 1 C * 6 C ' €■ 

n>l 

Thus we can easily compute a sequence bi, 62, ... of bounds such that Theorem [2] can be 
applied. To simplify notations let 

• ^mi,m 2 ,... - 

In the definition of above the coefficients of the corresponding power series in \E , °° (Jl°° ) 
can get larger than we considered so far because, to simplify things, we bound the real and 
imaginary parts of the n-th coefficient by m n , not the absolute value of the coefficient itself. 
Thus an additional factor v2 is sufficient in the following corollary. Due to the structure 
of our algorithm, however, we can get rid of this additional factor in implementations by 
restricting II 00 to a proper subclass by a simple test for the absolute value of the coefficients. 
We can summarize the essence of this section by the following corollary: 

Corollary 1. Let c > 1 and let, for n > 1, m n be some dyadic number such that 

n + 2 ^ n + 2 . n 
c ■ e • — ^— < m n < c ■ e ■ — h 2 n . 

Then 

(a) A = inf{A // |/GM/ 00 (n 00 )}, 

(b) we can compute on given dyadic complex number z € D, a E II 00 , n G N and r € 
Yfl(0;l) 

(i) f {n) {z) and 

(ii) upper bounds 



and 



11 / c • e c • e 

/i r := 



V 2 ' ' {l-rf + 2-(l-r)2 +2 



of sup{|/'(z)| I z G D r } and sup{|/"(z)| | z G D r }, respectively, 
up to any precision, where f '■= J ^?°°(ct). 



4. The A- value of a single Function 

In this section we will show how Aj can be approximated. More precisely, we will give 
lower bounds on Xf which are at the same time approximately upper bounds for A. This 
will be enough to compute A up to any precision in the end. To this end, let a\a2 ■ ■ ■ be a 
fixed sequence in IT 00 which will be fed to all algorithms considered in this section as input. 
Furthermore, we will denote the corresponding functions ^°°(ai . . .) and f \E , °°(ai . . .) by 
/' and /, respectively. 
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As in the case of Bloch's constant, it suffices to search large discs in the image of a 
proper sub-domain of the unit disc to approximate Xf. One advantage of this is that we have 
to search only on a bounded image for such discs. The following lemma can be proven by a 
simple transformation of holomorphic functions f(z) i— > ^f(r-z) (see e.g. |Conway 1978| ). 

Lemma 3. Let r with < r < 1 and a domain D with D r C D C be given. Then 

rX < Xf(D) < X f . 

Now, approximating the image f(O r ) can be done straightforwardly. Furthermore, with 
quite basic methods, the largest disc inside such an approximation could be found easily. 
We are doing exactly this by the e-covering grids which we introduce next. However, an 
approximation will not necessarily mean that we have an approximation of Xf : there indeed 
exist examples, where small changes of / can lead to a large change in Xf. That means that 



the homotopic methods used in |Rettinger 2008 cannot be used here. Instead, we will show 



that the approximations we get by e-covering grids for small e are still suitable to bound A 
from above. 

Definition 1. Let e > and a bounded subset A C C be given. Then an e-covering grid 
of A is a tuple (e, 5, G) where 0<<5<e/4is some dyadic number and G is a non-empty, 
finite subset of bTL x 5Z such that 

(a) A n (<JZ x 6Z) C G and 

(b) VzG G.3a e A.\z — a\ < e/4. 



Notice, that for any e > e' > and any e'-covering grid (e' ; 5, G) of a set A, the tuple 
(e, 5, G) is an e-covering grid of A. In this sense any e'-covering grid of a set A is also an 
e-covering grid. 

Following the notation of Section [2j 1(A) denotes the radius of the largest disc inside a 
domain A C C. Furthermore let 

lie, 5, G) := 5 + max min U — y| 

2SG ye(<5Zx<5Z)\G 

for an e-covering grid (e,6,G). Then the easy to compute value l(e,S,G) gives us an 
approximation of the largest disc inside the "covered set" 

D(e,5,G) := [j Bs E (z). 

More precisely we get the following result: 

Lemma 4. Let (e, 5, G) be an e-covering grid of A C C Then 

1(A) < l(e,5,G) <l(D(e,S,G)). 

Proof. We start by proving the left inequality: Let H) r (z) be some disc with D r (z) C A. We 
can assume that r > 5 because otherwise we have r < 5 < l(e,6,G) anyway. There exists 
some y € G such that \y — z\ < 5. Let furthermore x S bTL x b7L be some point with x $lG. 
Then we have ^\ x - y \(y) 2 ^ and thus D| x _ y | +< 5(z) <2 A. Therefore |x — y\ + 5 > r and the 
statement follows. 

The right inequality follows immediately from the following claim: 

Claim 1. Let x 6 G be given and r := 5 + mm z€ ^ZxSZ)\G \x — z\. Then D r (x) C D(e, 5, G). 
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Let us assume that there exist y D(e,5,G) = UzeG^eO 2 ) sucn that \x — y| < r. 
Then there also exists some point z' G 5Z x <5Z such that |x — z'| < r and |z' — y| < \/2<5. 
As \/2b < 3/4 • e, z' cannot belong to G because otherwise y would belong to U^ec ^e^)- 
Furthermore, as z 1 ^ x and x, z' G 57* x 57 there exists some z £ 57 x 57 such that 




Figure 1: Proof of Claim [T] 



\z — z'\ < \[2 ■ 5 and \z — x\ < \z' — x\ — 5 < r — 5, which means z G G by the definition of 

r and furthermore _ 

/2 3 
\z -y\<2- y/2- 5 < 2 • ^ ' ~ 



4 4 



□ 



in contradiction to y g" (J^eG 

Besides, an e-covering grid of /(O r ) can be easily computed: 

Lemma 5. For given dyadic numbers e > and < r < 1 as input, we can compute an 
e-covering grid of f(D r ). 

Proof. Following the notation of Theorem 2 we can compute a dyadic upper bound fj,' r on 
the maximum of the values with z G D r . Furthermore let 5q > be a dyadic lower 

bound on e/(16 • // r ) and 

G D :=D r n(fcZxW). 
Notice that OgG^i and thus Gd 7^ 0. Furthermore we can compute for every z G Gd some 
approximation <i z with |d 2 — f(z)\ < e/16. 

Then (e, e/4, G) is an e-covering grid of /(B r ), where 

U {y^j^xj^\\d z -y\< 3 -- £ -}. 



zeG D 



To see item (a) of Definition [l] let y G /(B r ) n (|Z x |Z) be given. Furthermore let 
x G D r such that /(x) = y. Then there exists some z G Gd with \z — x| < y/2 ■ 5d and we 
get 

\y-d z \ = \f(x)-d z \ 

< \f( x )-f( z )\ + \f(z)-d z \ 

< fjf r -\x-z\ + ^ 

< V2-^ + ^ 

<- 3 e 
— 2 ' 8 

To see item (b) of Definition [T] let now y G G be given. Then there exist z G Gd such 
that \d z — y\ < (3/2) • (e/8) and we have 

inf \y-a\<\y- f(z)\ < \d z - y\ + \f(z) -^|<|.£ + ± = £. 
oe/(B r ) 2 o lb 4 



□ 



ON COMPUTABLE APPROXIMATIONS OF LANDAU'S CONSTANT 



9 



Up to this point we have followed the naive way of computing large discs in approxi- 
mations of /(D r ). We have already argued that this will not necessarily give us appropriate 
bounds on Xf. The main step we will take next is to show that for suitable r > r and suit- 
ably small e we get that for any e-covering grid (e, 5, G) of /(D r ) we have D(e, 5, G) C /(Df ), 
which in the end allows us to show that the bound given by l(e,S,G) is actually not too 
bad. 

Lemma 6. On given dyadic number < r < 1 we can compute dyadic numbers e > and 
f with r < f < 1 such that for any e-covering grid (e, 5, G) of /(D r ) we have 

(J Bs e (z) C /(D r ). 

Proof. As / is not constant, we can compute p > and f such that < r < f < 1 and 
> P for all z G C with \z\ = r. Furthermore, following the notation of Corollary 1, 
we can compute a dyadic upper bound p!' f on the maximum of for zeDf. Let A be 

a dyadic number with 

< A < p/(4p") and 2 • A < f - r. 
Then we have |/'(^)| > p/2 for all z£lf\ Df_2A and, with r := r — A we have < r < 
r < r < 1. By choice of A we do have actually more: / is injective on Da/2( z ) f° r each 
z £ C with 

r- A/2 < |z| < r + A/2. 
Thus we can apply Koebe's 1/4-Lemma to get the following claim: 

Claim 2. For each z £ C with f — A/2 < |z| < r + A/2 and each y with |/(^) — y\ < p-A/16 
we have y G /(Dr)- 

Assume now, that there exists some zo with | zq \ < r — A/2 and some yo with | yo — 
f(zo)\ < p-A/16 and yo ^ /(Dr)- Then for the holomorphic function h : D — > C, determined 
by /i(z) := yo~ f{z) for all z G D, we have /i(z) 7^ for all z G Dr, i.e. takes its minimum 
on the boundary dDr = { z \ \ z \ = r}. Therefore there exists some z with \z\ = r and 

\h(z)\ < \y -f(z )\ <P-^- 
Thus by Claim [2] we have yo G /(D r ), i.e. we can extend Claim [2] as follows: 

Claim 3. For each z £ C with \z\ < r + A/2 and each y with \ f(z) — y\ < p-A/16 we have 
yG/(D r ). 

Let now e := p ■ A/16, (e, 6, G) be some e-covering grid of /(D r ) and x G (J 2gG D3 £ (z). 
Then there exists some z G G with \z — x\ < (3/4) • e and, by condition (b) in Definition 1, 
there exist y G /(D r ) such that \z — y\ < (l/4)-e. Thus we have \x — y\ < (l/4)-e + (3/4)-e = 
p ■ A/16 and, by Claim[3]we have x G /(Dr) which proves the lemma. □ 

Finally, we can combine the results of this section as follows: 

Corollary 2. Given «i«2 ■ • • in U°° and n G N, we can compute a dyadic number / such 
that (1 - 2~ n )X <l<X f , where / = / ^°°(Qia 2 . . .). 

Proof. Given some n G N we can compute by Lemma [6] first some e > and f < 1 with f > 
1 — 2~ n such that for any e-covering grid (e, 5, G) of /(D 1 _ 2 -n) we have D(e, 5, G) C /(D r ). 
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Then, by Lemma [5] we can compute an e-covering grid (e,5,G) of /(B]_ 2 -n) and compute 
I := l(s,5,G). Finally, by Lemma [4] and Lemma [3] we get 

(1 - 2~ n ) ■ X < A / (B 1 _ 2 -„) < I < l(D(e, 5, G)) < A/(B f ) < A/. □ 

5. The Main Theorem 

The proof of our main theorem can now be simply done by covering the space of A-bounding 
functions by neighbourhoods given by the algorithm of the previous section. 

Theorem 3. We can compute approximations to the Landau constant up to any precision. 

Proof. Let n 6 N be given. We proceed in steps t = 0, 1, . . . as follows. For each t we 
consider the set 

W t := {wl°° | w £ LI*}, 
i.e. the set of infinite words starting with an arbitrary (finite) word followed by infinitely 
many 1. Now we apply the algorithm given in Corollary 2 to each word u) in Wt and receive 
an value l(u>) such that 

(1 - 2" n )A < < Xj^y 

If for any oj 6 Wt information on any symbol outside the leading t symbols of u are asked 
by the algorithm, then we go to the next step t + 1. Otherwise / := inf wg vy t is an 
appropriate approximation of A, i.e. (1 — 2~ n )A < I < A. 

The latter is obvious, because, as none of the algorithms asks any information outside 
the leading words w 6 FT*, the results on wl°° equals the results on wv for all v £ IT 00 . 
Thus / is indeed the infimum over all approximations to A j j of all / G ^f°° (H°° ) . 

It remains, therefore, to prove that for each n there does indeed exist some t such 
that our algorithm stops after this step. This can be seen as follows: As the algorithm of 
Corollary 2 stops on all a G LI 00 , there exist for each a some n(a) such that no information 
is asked by the algorithm on any symbol beside the leading n(a) symbols of a. Thus the 
computation of the algorithm is identical for all words which coincides with a on the first 
n(a) symbols. This means that for each a there exists an open neighbourhood of a in 
LI 00 such that the algorithm is identical on any word in this neighbourhood. As IT 00 is 
compact, there is a finite sub-covering of LI 00 by such neighbourhoods. Let a ,...,a m be the 
corresponding words. Then the algorithm will stop in step 

t := max{n(a l ) | 1 < i < m} 

or before. □ 

The most interesting open problem is whether the conjectures on the Landau and Bloch 
constant hold. If so, the constant can clearly be computed in polynomial time. To this end, 
our algorithm and the algorithm given in [Rettingcr 2008] can present holomorphic functions 
which are very near the optimum for /3- and A-values, thus giving possibly new insights on 
the kind of functions involved. 

Concerning our algorithm the main intriguing problem is to improve the complexity 
bound, which is roughly double exponential, to an acceptable running time. The problem 
on reducing the time complexity of our algorithm is that the functions we have to consider 
can explode when reaching the boundary <9D, where evaluation can be quite expensive. 
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